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I review the problem of parton saturation and its implications through three in principal 
different approaches, but somewhat related: saturation in a geometrical approach, QCD 
saturation through the Color Glass Condensate and shadowing in a Pomeron approach. 



1. Introduction 

In the recent experiments like DIS at HERA or the heavy-ion experiments at RHIC, 
and also in expected LHC at CERN, the number of involved partons is very large, due to 
the high energy and/or the high number of participants of those experiments. These high 
parton densities should in principal lead to an extremely huge multiparticle production, 
but experimentally we have seen that this is not the case. So there should be a mechanism 
that reduces the number of created particle. Two kind of phenomena have been proposed. 
On one side there is the possibility of saturation in the initial state of the collision. On 
the othe side, there is the proposal of the creation of a Quark Gluon Plasma, and it 
has been presented as a final state interaction mechanism. Here, I review the problem 
of parton saturation and its implications through three in principal different approaches, 
but somewhat related: saturation in a geometrical approach, QCD saturation through 
the Color Glass Condensate and shadowing in a Pomeron approach. 



2. Geometrical approach to saturation 

2.1. String models and percolation 

In many models of hadronic and nuclear collisions, color strings are exchanged between 
the projectile and the target. Those strings act as color sources of particles through the 
creation of q — q pairs from the sea. The number of strings grows with the energy and 
with the number of nucleons of the participant nuclei. 

In impact parameter space these strings are seen as circles inside the total collision 
area. When the density of strings becomes high the string color fields begin to overlap and 
eventually individual strings may fuse, forming a new string -cluster- which has a higher 
color charge at its ends, corresponding to the summation of the color charges located at 
the ends of the original strings. The new string clusters break into hadrons according to 
their higher color. As a result, there is a reduction of the total multiplicity. Also, as the 
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Figure 1. From left to righ: Density of strings in the transverse space, from low energy 
and/or low number of participants to high energies and/or high number of participants. 
In the last circle we show percolation. 



energy-momenta of the original strings are summed to obtain the energy-momentum of 
the resulting cluster, the mean transverse momentum of the particles created by those 
clusters is increased compared to the one of the particles created from individual sources. 

As the number of strings increases, more strings overlap. Some years ago, it has been 
proposed in Ref. [ that above a critical density of strings percolation occurs, so that 
paths of overlapping circles are formed through the whole collision it is represented 

in Fig. ^ Along these paths the medium behaves like a color conductor. Also in [Q], we 
have made the remark that several fused strings can be considered as a domain of a non 
thermalized Quark Gluon Plasma. The percolation gives rise to the formation of a non 
thermalized Quark Gluon Plasma on a nuclear scale. 

Note that here we are not speaking about a final state interaction phenomenon, since 
there is no thermalization involved. In fact, what we are trying is to determine under 
what conditions the initial state configurations can lead to color connection, and more 
specifically, if variations of the initial state can lead to a transition from disconnected 
to connected color clusters. The results of such a study of the pre-equilibrium state in 
nuclear collisions do not depend on the subsequent evolution and thus in particular not 
require any kind of thermalization. 

The structural problem underlying the transition from disconnected to connected sys- 
tems of many components is a very general one, ranging from clustering in spin systems to 
the formation of galaxies. The formalism is given by percolation theory, which describes 
geometric critical behavior. 

2.2. Percolation theory 

Consider placing N small circular discs (color sources, strings or partons) of radius r 
onto a large circular manifold (the transverse nuclear plane) of radius R; the small discs 
may overlap. With increasing density 

A%r 2 . , 

" = w « 

this overlap will lead to more and larger connected clusters. The most interesting feature 
of this phenomenon is that the average cluster size increases very suddenly from very small 
to very large values. This suggests some kind of geometric critical behavior. In fact, the 
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cluster size diverges at a critical threshold value n c of the density. This appearance of an 
infinite cluster at 77 = n c is defined as percolation: the size of the cluster reaches the size 
of the system. r/ c has been computed using Monte Carlo simulation, direct connectedness 
expansion and other different methods. All the results are in the range n c = 1.12 1.175. 

In our model [ 1. we had proposed a fixed radius for the independent color sources 
of r = 0.2 4- 0.25 fm, that corresponds to a momentum around 1 GeV. This value has 
been obtained from Monte Carlo simulations in the framework of the String Fusion Model 
Code (SFMC) [|2] made at SPS energies. According to eq. (JTJ), in order to estimate the 
density rj, one needs to know the number of sources N. In our model, it is obtained 
from the SFMC, that, for nucleus-nucleus collisions, takes into account two contributions: 
one proportional to the number of participant nucleons -valence-like contribution- and 
another one proportional to the number of inelastic nucleon-nucleon collisions. Note that 
N will depend on the energy y/s (or equivalently, on x) and on the number of participant 
nucleons A, so in some way the condition to achieve percolation depends on A and s, 
rj = r](A, x). iiR 2 corresponds simply to the nuclear overlap area, Sa, at the given impact 
parameter. This overlap area can be determined in a Glauber study, using Woods-Saxon 
nuclear profiles. That leads to the following results: In our model, at SPS energies, the 
critical threshold for percolation could eventually been achieved for the most central Pb- 
Pb collisions, and for sure in Au-Au central collisions at RHIC energies and even in p-p 
collisions at LHC energies. 



2.3. The size of the color sources 

Let us now introduce a new question: If you look at a fast nucleon coming at you, what 
do you see? It depends on who is looking. Another nucleon sees a disc of radius r ~ 1 
fm and a certain greyness. A hard photon, with a resolution scale Q^ 1 << 1 fm, sees 
a swarm of partons. How many there are depends on the resolution scale: given a finer 
scale, you can see smaller partons, and there are more the harder you look. The partons 
in a nucleon have a transverse size tt determined by their transverse momentum kx , with 
r T ~ 1/Jzt ■ The scale Q^ 1 specifies the minimum k^ 1 resolved, so the probing photon 
sees all partons in the range < kx < Q, or equivalently it sees all the partons with a 
radius tt > 

So the partonic size is through the uncertainty relation determined by its average trans- 
verse momentum, r 2 ~ 1/ <k\>, for a given resolution scale, < k\ >~ Q 2 . 
These ideas are illustrated in Figs. El and El 

In order to know if the percolation density is achieved in a collision, we need to compute 
our eq. Then it is necessary to know the number of initially created sources (partons 
or strings, it depends on the kind of model we are using) N, and the size of those sources 
r. Remember that in the previous section we have used a fixed size of r = 0.2 4- 0.25 fm 
that corresponds to a moment of 1 GeV. In [Q] N is calculated in the framework of the 
SFMC. 

In [ 0] , the number of sources N is calculated through the Wonded Nucleon Model 
for nucleus-nucleus collisions, using P.D.F.'s for the nucleonic density. The authors use 
here an slightly different strategy: Instead of fixing the radius of the initially created 
sources, they estimate the momentum of those sources that will lead to percolation. They 
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Figure 2. From left to righ: The structure of an incoming nucleon seen for increasing 
resolution. 




Figure 3. The resolution scale. 
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transform then eq. into the following one: 
Nirr 2 N 



lc nR 2 QIB? 

The condition for percolation, taking r\ c = 1.12 will be then: 



(2) 



Q " ~ 1.12R 2 ■ (3) 

They find that for central Pb-Pb collisions at SPS energies, the critical momentum for 
percolation is Q 2 C ~ 1 GeV 2 , in accordance with our result. For central Au-Au collisions 
at RHIC, Q 2 ~ 2.5 GeV 2 . Note that again the condition to achieve percolation depends 
on A and \fs, Q c = Q c (A,x). 

Beyond the percolation point, one has a condensate, containing interacting and hence 
color-connected sources of all scales kr < Q. The percolation point thus specifies the 
onset of color deconfinement; it says nothing about any subsequent thermalization. 

3. The Color Glass Condensate 

Now we arrive to another approach, the QCD saturation through the formation of a 
Color Glass Condensate (CGC) [Ej. The idea is the following: At high energy, the QCD 
cross-sections are controlled by small longitudinal momentum gluons in the hadron wave 
function, whose density grows rapidly with increasing energy or decreasing x, due to 
the enhancement of radiative process. If one applies perturbation theory to this regime, 
one finds that, by resumming dominant radiative corrections at high energy, the BFKL 
equation leads to a gluon density that grows like a power of s and in consequence to 
a cross-section that violates the Froissart bound. Nevertheless, the use of perturbation 
theory to high-energy problems is not obvious. In fact, the BFKL and DGLAP equations 
are linear equations that neglet the interaction among the gluons. With increasing energy, 
recombination effects -that are non-linear- favored by the high density of gluons should 
become more important and lead to an eventual saturation of parton densities. 

These effects become important when the interaction probability for the gluons becomes 
of order one. Taking as the transverse size of the gluon and xG ^?p as the density 
of gluons, the interaction probability is expressed by 

<x a Nc xG(x,Q 2 ) 

Q 2 X nR 2 ' { ' 

Equivalently, for a given energy, saturation occurs for those gluons having a sufficiently 
large transverse size r\ ~ 1/Q 2 , larger than some critical value l/Q s (x,A). So the 
phenomenon of saturation introduces a characteristic momentum scale, the saturation 
momentum Q s (x,A), which is a measure of the density of the saturated gluons, and 
grows rapidly with 1/x and A (the atomic number). The probability of interaction -that 
can be understood as "overlapping" of the gluons in the transverse space- becomes of 
order one for those gluons with momenta Q 2 ^ Q s (x, A) where 

a xG(x,Q 2 ) _ (color charge) 2 

Q s {x, A) = a s N c — — = . (5) 

7i R z area 
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For Q 2 ^ Q 2 (x, A), the non- linear effects are essential, since they are expected to soften 
the growth of the gluon distribution with r = ln(l/x). For a nucleus, xG^(x, Q 2 ) oc A 
and tiR\ oc A 2 / 3 , so eq. (jSJ) predicts Q 2 oc A 1 / 3 . One can estimate the saturation scale by 
inserting the BFKL approximation into eq. (0). This gives (with 5^1/3 and A ~ ca s in 
a first approximation): 

Q 2 s (x,A) ~ A s x~\ (6) 

which indicates that an efficient way to create a high-density environment is to combine 
large nuclei with moderately small values of x, as it is done at RHIC. In fact the estimated 
momentum for saturation at RHIC will be Q s = 1 4- 2 GeV, in accordance with the result 
of the previous section. 

We call the high density gluonic matter at small- a; described by this effective theory 
a Color Glass Condensate: Color since gluons carry color under SU(N C ); Glass since 
we have classical coherent fields which are frozen over the typical time scales for high- 
energy scattering, but randomly changing over larger time scales. So we have a random 
distribution of time-independent color charges which is averaged over in the calculation 
of physical observables, in order to have a gauge independent formulation -in analogy to 
spin glasses-; and Condensate because at saturation the gluon density is of order l/a s , 
typical of condensates, so we have a system of saturate gluons that is a Bose condensate. 



4. Shadowing 

Note that in the string models, like Quark Gluon String Model or Dual Parton Model 
each string corresponds to the exchange of pomerons, and the interaction among the 
strings would correspond to interaction among the pomerons through the triple pomeron 
vertex. In particular, in the DPM, some kind of saturation is also included through the 
shadowing in the initial state of the collision [Ej. In fact, in absence of shadowing, the 
DPM at high energies leads to multiplicities that scale with the number of binary collisions 
rather than to a scaling with the number of participants. This is a general property of 
Gribov's Reggeon Field Theory which is known as AGK cancellation, analogous to the 
factorization theorem in perturbative QCD and valid for soft collisions in the absence of 
triple Pomeron diagrams. Because of this, a dynamical, non-linear shadowing has been 
included in the DPM. It is determined in terms of the diffractive cross-section. It is 
controlled by triple pomeron diagrams and it should lead to saturation as s — > oo. In 
fact, the shadowing corrections in this approach have a positive contribution to diffraction 
and a negative one to the total cross-section. 

The shadowing is controlled by the following quantity: 



rM 2 

f(x, Q 2 ) = 4vr / dM 2 
Jm 2 . 

mm 



1 da rv 



^y^nucleon dAd dt 



V 

F\{t min ) (7) 

i=0 



which is equivalent to (see below for details): 

da ppp 



rVrnax j 

F(s,y*)=^ dy — 

J Vmin 0~p\. 



dydt 



F^(t m i n ) } (8) 

i=0 
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where 

F A (tmin) = I d 2 b J (bV^t~)T A (b), (9) 



with t m in = — m%x'p and itln the nucleon mass. It represents the coherence effects, and 
at RHIC or higher energies can be taken equal to 1. 

Quantities as d ^ t \t=o or ^ ^^\ t=0 have been computed through F 2 (x, Q 2 ) 

and F$(Q 2 ,x P ,(3) data on DIS scatterring. 

The above expressions correspond to the case with only two scatterings. In order to 
include higher order rescatterings it is neccessary to do an unitarization. Two possible 
models apply here, eikonal unitarization or Schwimmer unitarization. Eikonal unitariza- 
tion will produce larger shadowing than Schwimmer, as can be expected by comparing 
the second non-trivial order in the expansion of both expressions. Both approaches have 
been checked with experimental data from DIS and heavy-ion collisions by Capella et al. 
and Frackfurt et at. Experimental data favorize the Schwimmer unitarization, that leads 
to: 

RiA/^eon = 1 + Afix ] Q2)TA{b) ■ (10) 

Note that shadowing in nuclei is usually studied through the ratios of cross sections per 
nucleon for different nuclei. This is the meaning of R(A/nucleon) above. 

Then the factor for reduction of multiplicities at fixed impact parameter b is: 

-lab{o) 

Ra{b)(P) is given by the r.h.s. of (jTUj) multiplied by Ta(b)(P) -nuclear profile functions-, 
2~ab(&) = J d?s Ta(s)T b (Ij — s), and with f(x,Q 2 ) expressed in terms of rapidity. 

As I said above, the shadowing is controlled by the quantity f(x,Q 2 ) (the shadowing 
grows with /). Their integrations limits are defined as: M^ in = 4m 2 = 0.08 GeV 2 , while 
the upper one is taken from the condition: 

XP = X I -Q2 J ^ Xpmax M max = Q 1 J, ( 12 j 

with x Pmax = 0.1 0.2. 

In order to expressed it as a function of the rapidity of the produced particles F(s, y*), 
we have two possibilities: the first one is using (J7J), but with the integration limits inspired 
by the parton model for hard processes: for projectile A (target B), 

x A (B) = ^fe ±r , (13) 
with y* > for the projectile hemisphere and y* < for the target one, and my = 



m 2 + j)\ the transverse mass of the produced particle. On the other hand, we can also 
compute the reduction factor from the formula 

rymax 1 da ppp 
F(s,y*)=A7T dy 



<7p(s) dydt 



Fi{t min ) } (14) 

t=o 
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where <Tp(s) is the single Pomeron exchange cross section and a dydt the triple Pomeron 
cross section. 

Using the standard triple Pomeron formula for the latter, we get 
1 da ppp 

= CAexp(Ay), (15) 



(Jp(s) dydt 



t=o 



with C = g pp(°^ pp (°) ; gpJO) the Pomeron-proton coupling and rppp(0) the triple Pomeron 
coupling, both evaluated at t — 0. That leads to the following values: C = 0.31 fm 2 , 
A = 0.13, fixed from DIS data. The same integration limits used above correspond to: 

(MB)) j ( s \ / R A m N \ {A (B)) = I ln ( J_\ * (m) 

where y* is the center of mass rapidity and m T the transverse mass of the produced 
particle. Finally we get: 

dy A exp (Ay) = C[exp (Ay max ) - exp (Ay min )\ . (17) 

5. Conclusions 

We have compared different models that takes into account saturation in different ways: 
from the semi-phenomenological fusing color sources picture for the soft domain including 
percolation, the QCD saturation through the Color Glass Condensate and those which 
follow from the pomeron approach and that takes into account shadowing in the initial 
conditions. 

In fact, it seems that the exchanged of elemental objects, color sources -strings, partons 
or pomerons-, should lead to a saturation in the initial conditions when the densities are 
high enough. In particular, in which concerns multiplicity dependence on the number 
of participants A, the cited approaches coincide in the results: when going to very high 
energies, the yield of created particles will be proportional to the number of participant 
nucleons when the shadowing/saturation effects are taken into account. In abscence of 
initial effect those yields behave as the number of collisions, so as A 4 / 3 . 
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